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Abstract 

Quantum string, which was brought into discussion recently as a model for the 
stripe phase in doped cuprates, is simulated by means of the density-matrix- 
renormalization-group method. String collides with adjacent neighbors, as it won- 
ders, owing to quantum zero-point fluctuations. The energy cost due to the collisions 
is our main concern. Embedding a quantum string between rigid walls with separa- 
tion d, we found that for sufficiently large d, collision-induced energy cost obeys the 
formula ~ exp(— ^d") with a = 0.808(1), and string's mean fluctuation width grows 
logarithmically ~ logd. Those results are not understood in terms of conventional 
picture that the string is 'disordered,' and only the short-wave- length fluctuations 
contribute to collisions. Rather, our results support a recent proposal that owing 
to collisions, short-wave-length fluctuations are suppressed, but instead, long-wave- 
length fluctuations become significant. This mechanism would be responsible for 
stabilizing the stripe phase. 

1. Introduction 

Recently, Zaanen brought up a problem of 'quantum string'0, which is a line-hke object 
subjected to line tension, and it wonders owing to quantum-mechanical zero-point fluctu- 
ations. Central concern is to estimate the interaction between adjacent strings as it won- 
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ders quantum-mechanically and undergoes entropy-reducing 'collisions' with their neighbors. 
Statistical mechanics of quantum-string gas would be relevant to the low-energy physics of 
the so-called stripe phase observed experimentally in doped cuprates&fl. In particular, one 
is motivated to gain insights how the stripe pattern formed in cuprates acquires stability. 
Actually, a good deal of theoretical analyses had predicted tendency toward stripe-pattern 
formation's. However, first-principle simulations on the t-J model still remain controver- 
sial about that issu The aim of the aforementioned recent analysis^ is to shed light, 
particularly, on the role of the quantum fluctuations and the entropy-reducing collisions at 
the expense of disregarding microscopic ingredients of quantum string. 

In the path-integral space-time picture, quantum string spans a world sheet as time 
evolvesEl. Hence, one may wonder that physics of quantum string might bear resemblance 
to that of membrane, and so quantum string is readily understood through resorting to the 
past outcomes about membrane&B. However, nature of elasticities is different: Quantum- 
string's elastic energy is quadratic in its slope (line tension), whereas membrane's is quadratic 
in its curvature (bending elasticity). This seemingly slight difference causes, according to 
Ref. 0, quite distinctive behaviors. 

When quantum string is laid down in a free space, its Hamiltonian is quadratic, and thus 
no peculiarity emerges. (We will introduce the Hamiltonian afterwards.) Yet, when it is 
squeezed by adjacent neighbors or walls, there immediately arise awful complications due to 
many-body correlations. In order to tackle with the string gas, Zaanen employed Helfrich's 
techniqueii'0, that has been utilized in the course of the studies of stacked membranes. 
Thereby, he arrived at the conclusion that collision-induced energy cost / would be described 
by / ~ exp(— Ac?^/^) with string interval d. That result cannot be understood in terms of 
conventional argumentH which gives / ~ exp(— Ad^); this argument is based on such picture 
that string is 'disordered,' and the collisions are mainly due to short-wave-length meandering 
modes. (We will outline this argument in Sec. |3 B| .) On the contrary, Zaanen emphasizes 
the role of long-wave-length fluctuations, because his theory is sensitive to infrared cutoff. 

The purpose of this paper is to judge the validity of those scenarios by performing first- 



principle simulations. We put a quantum string between rigid walls with spacing and 
measured its repelling interaction by observing pressure against the walls. (This trick has 
been used in the studies of the fluctuation pressure of (classical) membraneEHli.) The 
Hamiltonian, which we had simulated, is given by, 

^ = E + nx^)) + E |(^^ - ^m)'. (1) 

Here, Xi denotes the operator of transverse displacement of a particle at ith site, and pi is 
its conjugate momentum. They satisfy the canonical commutation relations = ih6ij, 

[xi,Xj] = and [pi,Pj] = 0. V{x) is rigid-wall potential with spacing d; 

for < x < 

V{x) = l - - . (2) 

I oo otherwise 

S denotes line tension which puts particles into line. Classical version of this Hamiltonian has 
been used as a model for line dislocations and steps on (vicinal) surfaceJ^HH. Note that for 
sufficiently large S, one can take continuum limit, with which one arrives at field-theoretical 
version of quantum string. Such continuum- limit version was studied analytically in Ref. |l|. 
For small S, however, continuum-limit treatment fails, and our lattice- version Hamiltonian 
exhibits characteristics which are not considered in the previous study. 

The rest of this paper is organized as follows. In the next section, we explicate our 
simulation scheme. We used diagonalization method: Note that elastic models such as ours 
(0) have vast number of vibration modes, which overwhelm computer-memory size. Empha- 
sis is laid upon the point how we had adopted the idea of density-matrix renormalization 
groupSii to quantum string so as to reduce the number of Hilbert-space bases. Our algo- 
rithm is indebted to recent develop mentsS, where a polaron model (lattice vibration) was 
simulated with use of the density-matrix renormalization group; see also Ref. Following 
the preparations, in Sec. ^, we perform numerical calculation. We show that our first- 
principle simulation supports aforementioned Zaanen's scenario. In the last section, we give 
summary and discussions. 
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2. Details of numerical method: density-matrix renormalization group 

In this section, we explain our simulation algorithm. Emphasis is laid on the point how 
we had adopted the idea of the density-matrix renormalization groupH'll in order to treat 
quantum string. Afterwards, we will show preliminary simulation data so as to demonstrate 
its performance. 

Quantum string is made of many particles connected with line tension S. Each particle 
spans infinite-dimensional Hilbert space. Hence, one is forced to truncate, somehow, the 
number of bases in order to diagonalize the Hamiltonian. Even though one truncated bases 
of each particle, the total bases of quantum string, as a whole, would exceed available 
computer-memory size. This difficulty arises inevitably in treating elastic (bosonic) degrees 
of freedom by means of diagonalization. Recently, however, Zhang et al. reported that 
the difficulty is overcome by an ingenious application of the density-matrix renormalization 
grou p§. Actually, they demonstrated that a polaron model is diagonalized successfully; see 
also Ref. In the following, based on these recent developments, we propose a scheme to 
diagonalize the quantum-string Hamiltonian (|l|). 

Our simulation algorithm proceeds recursively; after the completion of one recursion 
(renormalization), the number of treated particles (length of quantum string) is enlarged by 
two, and subsequent renormalization follows. In the following, we explain steps constituting 
one renormalization procedure: Suppose that our system (quantum string) is decomposed 
into four parts, that is, block, site, site and block, in this order. After one renormalization 
is completed, composite systems of block and site are 'renormalized' into a new block with 
block Hilbert-space dimension unchanged. 

As would be guessed, 'site' merely stands for one particle of quantum string. We then 
need to prepare bases to represent the Hilbert space of 'site.' We had chosen such M bases 
that are the eigenstates \n), (n = 1 ~ M) of the intra-site Hamiltonian Ti, = p^/(2m) + V(a;) 
with energy En = h^n'^n'^/ {2mcP'). (Therefore, the direct product of \n), gives the eigenstate 
of the total Hamiltonian 7i, provided that the line tension is tuned off (S = 0).) With respect 
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to these bases {l^i),} {n — 1-M). we represent the matrix of the Xi operator; 




,{n\x\m), = < 



TT- 



2d 



i {^^d^ - (;^^) ^) = (even, odd) or (odd, even) 



(3) 



nm 







otherwise 



The truncation bound M is one source of numerical errors. We need to choose M sufficiently 
large; we will demonstrate afterwards that this truncation does not deteriorate simulation 
precision in practice. 

Let us turn to explaining 'block.' The 'block' stands for a part (fragment) of quantum 
string, and actually, it contains many particles in it. Hence, at a glance, one may wonder that 
the Hilbert-space dimension of block would be extremely large. Yet, owing to the density- 
matrix renormalization, the dimension is reduced so that it can be stored in computer 
memory. The Hilbert space of block is spanned by those bases \n)B {n — 1 ^ m). The bases 
are to be prepared in the preceding renormalization procedure; see below. With respect to 
these bases, one has to represent the intra-block Hamiltonian He and the coordinate operator 
of the particle at the end of 'block' xb] see below. At the initial stage of renormahzation, 
the block is merely a 'site,' and so we start with Hb — Ti-,, xb — x, and m — M. 

Provided that the above matrices are at hand, the (total) Hamiltonian of quantum string 
is expressed in terms of them; 



ijkl 

With use of ipijki, we obtain the density matrix for the left-half sub-system (B -|- •); 
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+i®i®n,^i + i®i®i®nB 



+—{xl 0i(8)i0i-2xB®x. (8)i(8)i-Fi(g)a;^(g)i(g)i) 
-F-(i ® x\®\®\ - x,® x,®l^\®\® x\®\) 
+— (i ®1®x\®\-21®\®x,®xb + ^®^®^®x^. 



(4) 



Diagonalizing this matrix, one obtains the ground state. 



W) = II^ijfciN)B|j).|^).|OB- 



(5) 



[P\i,r.i'd' = Y.^iJki'^h'ki- (6) 

kl 

Diagonalizing this, we obtain the eigenstates and the eigenvectors; p\un) = Wn\un) with 
Wi > W2 > ■ ■ ■ > W(^Mxm)- Those bases \un) with large weight Wn would be important 
(relevant) in order to describe the physics of the sub-system of block and site. Therefore, we 
store the bases with n = 1-m, and discard the others. This criterion is the essence of 
the so-called density-matrix renormalization groupS^l. That truncation may cause another 
source of numerical errors. Later, we will demonstrate that this error is very small. Finally, 
we perform the 'density- matrix renormalization'; 



['HB']nm — \ Ur 



Kb ® i + i ^ n. + - (xl (g) i - 2xb ^ X, + i ® xl 
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Um ) (7) 



X, 



Um) ■ (8) 



Now, a renormalization is completed. We can restart the next renormalization from the 
beginning, replacing the renormalized block B' with B. It is to be noted that through the 
renormalization, the block dimension is kept within m. 

In the following, we will show that the above algorithm actually works. In Fig. |T], we 
plotted the relative error of the ground-state energy 6Eg/Eg for the system with S = 4 
and d = 4. (This parameter condition is of great physical significance as would be shown 
in Sec. ^.) The system size is L = 8, and the dimension of 'site' is M = 9. Therefore, 
the full number of bases is = 43046721, which is about to exceed the limit of available 
computer-memory size and is barely manageable with full diagonalization method. For 
these full bases, we calculated the 'exact' ground-state energy, while 'approximate' energy 
is calculated by means of the density-matrix renormalization group with truncated 'block' 
dimension m. The energy difference between them gives SEg. The relative error 6Eg/Eg is 
plotted in Fig. ^ We achieve very small error 10"*^ with m = 20, for which the total number 
of bases is no more than 20^ ■ 8^ = 25600. Hence, we see that our algorithm works efficiently. 
The inset of Fig. |I] shows the distribution of the density-matrix eigenvalues {wn}- Usually, 
Wn is utilized for monitoring 6Eg/Eg, because they look alike. However, in our case, there 
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are discrepancies in their magnitudes. The discrepancies may be, perhaps, due to the fact 
that our Hamiltonian-matrix elements distribute over wide range. By the way, we chose 
m corresponding to Wm = 10~^^ in Sec. ^. In this way, we kept precision within 10~^. 
Typically, we need, at most, m = 30 bases. 

In the above, we have checked that the truncation of block dimension m does not harm 
any reliability of our simulation. Finally, we will examine the influence of site-dimension 
truncation M. In order to elucidate that, it is sensible to monitor, 

= Y.^ijkl^ij'kl- (9) 
ikl 

As would be apparent from the definition, [p,]nn tells the degree of significance of the state 
|n),. Because we use the bases of n =1-M and discard the others, it should be checked 
whether [p,]mm is small enough. We see from Fig. 0, that the M = 9 state is of very rare 
probability 10~^. In the subsequent simulations in Sec. ^ we impose even severer request 
[p,]mm = 10^^°. In order to match this request, we need, at most, M = 20 bases; hence, 
the maximal total number of Hilbert-space bases is no more than ■ = 30^ • 20^ = 
360000. As is mentioned above, the full diagonalization of Figs. requires the Hilbert- 
space dimensions 43046721 for string length L = 8. It is far beyond the capability of 
the diagonalization method to treat longer string, unless we resort to the density-matrix 
renormalization group. 

3. Numerical results and discussions 

In this section, we present simulation results. In our simulation, we treated sufficiently large 
system sizes by repeating renormalizations, until simulation result converges to a (thermo- 
dynamic) limit. System parameters of m and h are fixed; namely, we set m = 1 and h = 1. 
those parameters just fix the coefficient of the kinetic-energy term. Therefore, the choice of 
parameters does not harm any generality. Technical details are to be referred to the previous 
section. 
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A. Collision-induced energy cost / and elasticity modulus B 



In Fig. ^ we plotted the collision- induced energy cost, 

y ^ E^{d) - Eg{d ^ oo) ^^^^ 

for the system with line tension S = 4 and various wall spacing d. Here, Eg{d) denotes 
the ground-state energy per one particle for wall spacing d. As would be apparent from 
the definition (|TUp, / measures excess energy cost (per unit volume) due to the presence of 
walls. It is notable that Eg{d — > oo) is calculated exactly, 



EJd ^ oo) = — / dk\ — 



k 

sm — 
2 



(11) 



because for ^ oo, the Hamiltonian reduces to quadratic form. 

In Fig. 1^, we notice that two distinctive regimes exist: For d < 2, the collision energy / 
decreases obeying power law, 

/ ~ l/d\ (12) 

whereas for d > 2, it decays 'rapidly'; actually, it is our main concern to clarify how 'rapid' it 
decays. The behavior for < 2 is understood immediately: Suppose that the line-tension is 
turned off (E = 0), each particle becomes independent, and it reduces to a text-book problem 
of 'particle in a box,' for which the ground-state energy is solved exactly; Eg = U^ti'^ / (2md'^). 
Hence, we arrive at the relation / = Eg/d ~ 1/d^. To summarize, for small d < 2, inter- 
particle interaction is irrelevant. On the other hand, for large d, the inter-particle interaction 
may become relevant. Because of this, the particles become correlated, and / drops very 
rapidly. 

In order to elucidate characteristics of fluctuation-induced interaction, it is sensible to 
calculate the elasticity modulus which is defined by the formula, 

B ^ (13) 
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In Fig. ^, we plotted B for the same parameter range as that of Fig. |^. As well as f, B also 
exhibits two distinctive regimes. For d < 2, we see that B decays in the form B ~ 
This behavior is understood immediately, if we remember / ~ l/rf^. The behavior for large 
d > 2 is still remaining mysterious. We will explore this in the next subsection. 



B. Scaling analyses 

Here, we re-examine the data of / and B presented in Figs. § and ^ For that purpose, 
we need to resort to past outcomes on / and B. First, we summarize recent remarkable 
predictions by ZaanenEI. In his theory, B works as a mean field, which is to be determined 
self-consistently. The self-consistency equation to be solved is. 

Appearance of logarithmic term brings about subtleties. It is hard to solve this self- 
consistency equation. However, asymptotic form for ^ oo is calculated as follows; 

/ - e--"'-'' (15) 
Br^d\-^"'''\ (16) 



In addition to Eq. ([1^) , under the assumption that the string- meandering modes are 'com- 
pactified' (frozen), he found another self-consistency condition. 



Note that the validity of this relation is checked, if we assume Eq. (12). Hence, this relation 
(l|) may be realized for small d. 

Besides those recent treatments, there exists an ingenious argument^! for tackling with 
entropic interaction. The argument, applied for quantum string (|I]), yields predictions dif- 
ferent from the above. Below, we outline this argument: As is mentioned in Introduction, in 
the path-integral picture, quantum string bears resemblance to (classical) membrane. When 
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quantum string is embedded in a free space, we can solve its mean deviation ~ (log/)^/^ with 
membrane's linear dimension /. Suppose that the wall width (mean membrane interval) is 
d, with use of this relation, one obtains an estimate of the surface 5* per one collision such 
as d ~ (log "s/S*)^/^. Assuming that each collision (contact-or-crossing) gives rise to entropy 
loss ~ /cBlog2, we obtain the collision-induced energy gain per unit surface; 

~ A;Blog2e-^"'^'. (18) 

The argument is based on the assumption that the string would be disordered spatially 
as in the Einstein-like view of a crystal. In other words, short-wave-length fluctuations 
contribute to collisions. On the contrary, Ref. emphasizes the signiflcance of long-wave- 
length fluctuations. Our aim is to judge which picture is valid by means of flrst-principle 
simulations. 

Guided by the above ideas, we will carry out scaling analyses. First, in Fig. ^, we plotted 
fd^/'^/^/B with use of the data shown in Figs. || and |[ Note that fd?/'^/^fB should be 
constant, if Eq. (|T^ holds. The scaling plot indicates that for < 2, in fact, the scaled data 
are kept constant. On the contrary, for d > 2, the scaled data drop suddenly, suggesting 
that Eq. ( P^D does not hold any more. In consequence, we conflrmed that for small d, 
string-meandering modes are 'compactifled' (irrelevant). These observations are consistent 
with those found in the previous subsection. 

Secondly, let us turn to large-rf regime. In Fig, ^ we plotted fd'^/B against d/B with use 
of the data shown in Figs, ^and ^. Note that the scale of abscissa is logarithmic. Therefore, 
the data should align, if the relation ([T^) is satisfled. In fact, for scaling regime d/B > 1, 
we see that the scaled data get aligned, and exhibit a positive slope. This result justifles 
the validity of Eq. (0). Namely, the mean-fleld treatment of Ref. ^ appears to capture the 
essence of this physics. For exceedingly large d, eventually, numerical data become scattered 
(unstable) gradually. This is a symptom of numerical errors. Note that for exceedingly large 
d, the elasticity modulus B is very small so that it suffers from tiny numerical errors. Hence, 
we cannot continue simulations for extremely large d. 
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In Fig. 0, we presented similar scaling data but for various S = 2, 4, 6 and 8. The plots 
show that the scaling-plot slopes are identical. This fact tells that the constant C in Eq. 



(|1^) is indeed universal with respect to S. It is, however, suggested that C is subjected to 
a correction to scaling, because plots do not overlap precisely. We had found that scaling 
data for very small S < 1 are not described by Eq. (p!^. As a matter of fact, the slopes 
are almost vanishing. Breakdown of Eq. (|14D for very small S is reasonable, because the 
equation is derived under continuum-limit treatment which is not justified for very small E. 
Therefore, such the region of small S lies out of the scope of Eq. (p!4|). By the way, we stress 
that our first-principle simulation does cover any conditions. 



C. Asymptotic form of / 

In order to confirm the above observation, we investigate the asymptotic form of / for 
d ^ oo. As is mentioned in the previous subsection, / should obey the asymptotic form, 

/ - e-^"^^ (19) 

with the exponent either a = 2/3 ( |T5|) or a = 2 ([T8|) . We expect that the former would 
be realized, because it is derived from Eq. (p!4|), whose validity is checked in Sec. |3 Bj We 
calculated a by means of the formula, 



I0g(l0g/(rfi)/l0g/(d2)) 

" = i — 7T7T\ ' y^^i 

with respect to adjacent two data points of d = di and d2 depicted in Fig. |. We plotted 
a in Fig. p. The scale of abscissa l/((logc/i -|- log(i2)/2)^ is chosen so as to achieve straight 
data alignment. We employed the least-square method in order to extrapolate the result for 
d oo, and thereby, we obtained the estimate a = 0.808(1). This result clearly supports 
a = 2/3 (|T3|) rather than a = 2 (|TB|). We notice that the convergence speed is rather slow; 



note that the abscissa scale of Fig. |^ is logarithmic. Hence, we found that the asymptotic 
form (|TB|) is realized for extremely large d. 
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D. Mean fluctuation width A 



So far, we have confirmed the vahdity of the relations (H) and (p!5[). Underlying physics 



of these relations are rather astonishingp: It is speculated that collisions rather contribute 
to straightening the quantum string. In this subsection, we will examine this remarkable 
scenario through observing the mean fluctuation deviation. 



A = v/(^fy^^. (21) 

In Fig. ^, we plotted A for the same parameter conditions as in Fig. ^. To our surprise, A 
grows logarithmically (A ~ log d) for large d. That is, A is much less than the wall spacing d. 
This feature contradicts our intuition that string would be short-range disordered, and the 
collisions would be due to microscopic fluctuations. Because Eq. (|^) relies on this picture, 
it turned out consequently, that the picture is incorrect for quantum string. Rather, our 
result supports the aforementioned speculation that string acquires macroscopic stiffness. 



4. Summary 

We have investigated fluctuation-induced repelling interaction of quantum string described 
by the Hamiltonian (|l]). First, we have developed simulation scheme based on the idea of the 
density-matrix renormalization group. We found that the scheme works very efficiently. As 
is demonstrated in Fig. |l|, ground-state energy is precise up to the eighth digit. Precision is 
crucial in our study, because we need to calculate second-order derivative in order to obtain 
the elasticity modulus B (p!3[). Secondly, based on those preparations, we have performed 
extensive simulations. Simulation data suggest that two distinctive regimes exist: For small 
d, intra-particle interaction dominates so that simulation data are understood by ignoring 
line tension (S = 0). That is, collision-induced energy cost / and the elasticity modulus 
B obey simple formulae such as /, -B ~ d~^. For large d, on the other hand, particles get 
correlated by line tension S, and physics becomes much harder to interpret. We made trials 
of several types of scalings in Figs. From those scaling plots, in consequence, we found 
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that the data of / and B agree with Zaanen's self-consistency condition ([M]). Moreover, we 
investigated the asymptotic form of /. We found / ~ exp(— C'd") with a = 0.808(1), which 
again supports Zaanen's result (|T5|) rather than Eq. (|18|). 

According to him, underlying physics of Eqs. (0) and (^) is quite peculiar; colli- 
sions with adjacent neighbors rather suppress short-range fluctuations, and consequently, 
quantum string is straightened macroscopically. In fact, we found that in Fig. P, string's 
fluctuation deviation A is bounded within A ~ log d, which is far less than the wall spacing 
d. The above observations tell that in essence, quantum-string's ground state is governed 
by infrared fluctuations. In this respect, it is decisively important to treat long quantum 
string; otherwise we could not have observed those features mentioned above. Hence, our 
simulation scheme is advantageous. 

We have conflrmed that owing to collisions, actually, order out of disorder sets in. This 
observation immediately leads an expectation that the stripe pattern observed in doped 
cupratesi^ is stabilized by this mechanism. To verify this deflnitely, one may need to 
explore 'stacked' strings. As for stacked membranes, it has been known that A^-membrane 
physics is essentially the same as that of a single membrane conflned between wa 
That is, there is a nice relation, with which one readily obtains A^-membrane behavior with 
use of smf/Ze-membrane result. It remains for future study to verify that similar relation 
holds as well as for stacked strings. 
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FIGURES 
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m 

Fig. 1. Relative error of the ground-state energy for the system with (i = 4, S = 4, L = 8 and 
M = 9. 'Exact' energy is calculated with respect to the full Hilbert-space dimensions = 9^, 
while 'approximate' energy is calculated for truncated 'block' bases m with use of the density-matrix 
renormalization group. Ground-state-energy difference between them gives bE^. Note that trun- 
cated-base calculation with small m reproduces full-diagonalization result very precisely. Inset 
shows the density-matrix eigenvalues {t^n}, which are used for monitoring simulation precision; see 
text for details. 
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Fig. 2. Probability weight [p,]nn of the state |n), is plotted. This plot indicates how 
site-dimension truncation M affects the reliability of simulation. Probability weight for n = 9 
is negligibly small; [p.Jgg ~ 10^^. Hence, truncation of those states with n > 9 does not deteriorate 
any reliability in practice. In Sec. ^ we impose even severer request \p^mm ~ lO^^'^. 
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d 

Fig. 3. Excess energy cost due to collisions / (|To| ) is plotted for S = 4 and various d. We see 
that for fi < 2, it obeys power law / ~ while for c? > 2, it drops rapidly. 
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Fig. 4. Elasticity modulus B (^) is plotted for the same parameter range as that of Fi^ 
We see that for d < 2, it obeys power law B ~ d^^, while for d > 2, B drops rapidly. 
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Fig. 5. jSl'^jB^I'^ is plotted with use of the data shown in Figs. |3| and ^ This plot shows 
that the relation ([T7|), which is valid for 'compactified' string, is realized for small d. On the other 
hand, for d > 2, it does not hold at all. 
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Fig. 6. fd^/B is plotted against d/ B with use of the data shown in Figs. ^ and ^. For d/ B > 1, 
the data approach to a straight line asymptotically. This result indicates that the relation ( [I^ ) is 
satisfied for the scaling region d/B > 1. 
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Fig. 7. Similar scaling analysis as Fig. ^ but for various S = 2-8. In the scaling regime 
dH/B > 3, slopes appear to be identical, suggesting that the constant C in Eq. (|l^) is universal 
with respect to arbitrary E. 
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Fig. 8. Exponent a (^) is plotted with use of the data shown in Fig. ^. With the least-square 
method, we extrapolated the data so as to obtain a = 0.808(1) for d oo. 
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Fig. 9. Meandering- fluctuation width (21) is plotted for the same parameter range as that of 
Fig. ^. For d > 2, A grows logarithmically (A ~ logfi). That is, string's fluctuation is far less 
than the wall spacing d. This result indicates that the string is straightened by collisions. 
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